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ABSTRACT 



Context. CoRoT and Kepler missions are now providing high-quality asteroseismic data for a large number of stars. Among 
intermediate-mass and massive stars, fast rotators are common objects. Taking the rotation effects into account is needed to cor- 
rectly understand, identify, and interpret the observed oscillation frequencies of these stars. A classical approach is to consider the 
rotation as a perturbation. 

Aims. In this paper, we focus on gravity modes, such as those occurring in y Doradus, slowly pulsating B (SPB), or Be stars. We aim 
to define the suitability of perturbative methods. 

Methods. With the two-dimensional oscillation program (TOP), we performed complete computations of gravity modes - including 
the Coriolis force, the centrifugal distortion, and compressible effects - in 2-D distorted polytropic models of stars. We started with 
the modes £ = i,n = 1-14, and t = 2-3, n = 1-5, 16-20 of a nonrotating star, and followed these modes by increasing the rotation 
rate up to 70% of the break-up rotation rate. We then derived perturbative coefficients and determined the domains of validity of the 
perturbative methods. 

Results. Second-order perturbative methods are suited to computing low-order, low-degree mode frequencies up to rotation speeds 
~100 kms"' for typical y Dor stars or ~150 kms"' for B stars. The domains of validity can be extended by a few tens of kms"' 
thanks to the third-order terms. For higher order modes, the domains of validity are noticeably reduced. Moreover, perturbative meth- 
ods are inefficient for modes with frequencies lower than the Coriolis frequency lil. We interpret this failure as a consequence of a 
modification in the shape of the resonant cavity that is not taken into account in the perturbative approach. 

Key words. Asteroseismology - Stars: oscillations, rotation - Methods: numerical 



1. Introduction 

CoRoT (C o nvecti on, Rotation and planetar y Tra nsits, 
iBa^lin et all '2006) and Kepler (iBorucki et alJ 12007 !) are 
space missions providing uninterrupted high-quality photome- 
try time series over several months or years ideally suited for 
asteroseismic study. Asteroseismology provides very accurate 
determinations of the stellar parameters (mass, radius, age, 
etc.) and probes stellar structure to constrain physical processes 
occurring in stars. The first step towards this goal requires 
correctly understanding the structure of the observed oscillation 
spectra, and especially coiTectly identifying the observed 
m odes. In the c ase of main-sequence (e.g. Michel et al. 20()a 



Here, we are concerned with gravity modes, i.e., low- 
frequency modes driven by the buoyancy force. They are excited 
and observed in a broad panel of stars, for instance, in y Dor, 
SPB and some Be stars. The y Dor stars form a class of main- 
sequence st ars with type aroun d FOV that can sometimes rotate 
rapidly (e.g. lDe Cat et al.ir2006l and reference therein), while the 
rotation rate of Be stars is extreme, usually very close to their 
break-up hmit n =i Q.k = ^JgmJ¥ (e.g. lFrematet al .1120061) . 

The effects of rotation on the oscillation modes can be treated 
as a perturbation where the rotation rate is the small parame- 
ter A 1 st-order correction has been proposed by Ledoux( 1951!), 
2nd-order by ISaiol (Il98ll) . iDziembowski & Goodd (11992). or 



Benomaretal. 2009; Ch aphnetalJ l2010l) and giaii t (e.^ 



tsenomar et at. zuuv; cn apunetatj l/UtU|) ana gian t (e.g . 
Migho et aL .2009; .Hekker et all 120091: iBedding et al.1 12010) 



FGK stars with solar-like oscillations, the spectrum structure is 
well understood, which eases interpretation. 

The spectra of classical pulsators is often noticeably more 
complex. For instance, the high-qual ity observations of 6 
Scuti (e . g. iGarciaH ernandez et al. 20091: iPoretti et alJl2009h and 
y Doradus (e.g. Mathi as et al.. .2009) stars have exhibited very 
rich and complex spectra of acoustic (p) and gravity (g) modes, 
respectively, containing several hundred -or more- modes. 
Interpreting their spectra is very challenging today. Indeed, these 
stars generally spin rapidly, so the effects of rotation on the mode 
frequencies must be considered. 



ISuarez et al.l (120061). and 3rd-order terms have been developed 
bv lSoufi et al.l (Il998l) . While perturbative methods are expected 
to be accurate enough for slowly rotating stars, their true domain 
of validity cannot be determined in the absence of exact calcula- 
tions to compare them with. 

In the past few years, calculations of p modes with both the 
centrifugal distortion and the Coriolis forc e have been perfo rmed 
in po lytropic models of stars dLignieres et al. 2006; Rees e et al. 
12006) and reahstic 2-D stellar structu res (Lovekin & D eupree 
200 ^Reese et al.ll2009h . lLignieres eTa l. (2006) and Rees e et al. 
(2006) have shown that, above Q ~ O.lSQ/f, perturbation meth- 
ods fail to reproduce low-degree and low-order p-mode frequen- 
cies (f < 3 and n < 10) with the accuracy of CoRoT long runs. 
The structure of the modes is also drastically modified, and this 
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leads to deep changes in the structure of the p -mode spectrum 
(iReese et al.l boOS' LJg nieres & G eorgeot'200g.'2009^. 

We used an oscillation code based on Reese et al. ( 2006) to 
perform g-mode calculations with a complete description of the 
rotational effects on the modes. In this paper, we focus on the 
limits of validity for perturbative methods. The models and the 
method are described in Sect. |2] We then derive the perturbative 
coefficients from the complete computations (Sect. [3]), and com- 
pare the results obtained with both methods to determine and dis- 
cuss the domains of validity for perturbative methods (Sect. |4]i 
before concluding in Sect.|5] 



2. Models and methods 

We consider fully radiative stars for this work. Since the grav- 
ity modes are driven by the buoyancy force, they cannot exist 
in convective regions. SPB and y Dor stars have large radiative 
zones with a convective core, and even a thin convective enve- 
lope for the latter. The effects of convective cores are not consid- 
ered here, since we are mainly interested in the general behavior 
of g modes under rotation effects. 

2A. 2-D stellar models 



As in iLignieres et al.l (2006) and iReese et al] (l2006l) . we ap- 
proximate the equilibrium structure of rotating stars with self- 
gravitating uniformly-rotating poly tropes. They are described in 
the co-rotating frame by the three following equations: 



Po 

Ai/To 



1 + 1//J 



Kpo 

Pogo 

AnGpo 



(1) 
(2) 
(3) 



where p„ is the pressure, p„ the density, 4'o the gravitational po- 
tential, K the polytropic constant, ji the polytropic index, G the 
gravitational constant, and g,, the effective gravity defined as 



go 



-^i^o 



■n^s^/2) 



with s the distance to the rotation axis. Due to the centrifu- 
gal distortion, the star is not spherical and a suited surface- 
fitting spheroid a l sys tem of coordinates (^, 6, (p) based on 
iBonazzola et al.l (Il998h has been used. Hereafter, we also clas- 
sically denote r the distance to the center, and z the coordi- 
nate along the rotation axis. This equation system is numerically 
solved with the ESTER (Evolut ion ST Ellaire en Rotation) code 
as described in .Rieutord et aP (12005!) . This is a spectral code 
using Chebychev polynomials in the ^-direction, and spherical 
harmonics Y^ with even { and m = in the horizontal one. We 
computed models decomposed on spherical harmonics up to de- 
gree Lmodei = 32. This resolution is high enough to accurately 
model the centrifugal effects for the maximal value of Q that we 
have considered. In the pseudo-radial direction, the resolution is 
the same as the one we use for the frequency computation (see 
Sect. 12.3! rir = 96 generally). 

To approximate a fully radiative star, we chose the polytropic 
index fi - 3. We considered models spinning with rotation fre- 
quency Q. between and QJD.K, where Q^: - ^GMjReq is the 
Keplerian break-up rotation rate for a star of mass M and equa- 
torial radius Rgq. 



2.2. Linearized equations for the oscillations 

In the co-rotating frame, the equations governing the temporal 
evolution of small adiabatic inviscid perturbations of the equi- 



librium structure read, in the co-rotating frame, 

d,p = -V-(p„v) 
p„d,v = -Vp + pgo - Po'^^ - 2pon X V 



d,p-cld,p = 



WgoW 



-V- go 



Aif/ = A-nGp 



(5) 
(6) 

(7) 

(8) 



where p, p, v, and i// are the Eulerian perturbations of density, 
pressure, velocity, and gravitational potential c^ = Txpolp,, the 
adiabatic sound speed and No the Brunt-Vaisala frequency, de- 
fined as 



Po 1 1 



Po 



(9) 



Fi = (51n/7/(91np)ad denotes the first adiabatic exponent. In 
Eq. (J?) we have used the structure barotropicity, ensured by the 
uniform rotation. 

The 2-D distribution of the Brunt-Vaisala frequency is 
shown in Fig. [1] for the most rapidly rotating model we have 
considered, together with the profiles of A^,, along the polar and 
equatorial radii, which are compared to the N„ profile of the non- 
rotating star Within about the inner half of the star, the devi- 
ations from sphericity induced by the centrifugal force remain 
limited. We also see that A^,, diverges at the surface of the poly- 
trope because po and po vanish there. 

Looking for time-harmonic solutions oc exp(/wO of the sys- 
tem ©-(O, we obtain an eigenvalue problem, which we then 
solve using the two-dimensional oscillation program (TOP). Th e 
details of this oscillation code closely follow lReese et al.l(l2006h . 
The equations are projected on the spherical harmonic basis Y^. 
Due to the axisymmetry of the system, the projected equations 
are decoupled relatively to the azimuthal order m, but in con- 
trast to the spherical non-rotating case, they are coupled for all 
degrees I of the same parity. 



(4) 2.3. Resolutions and method accuracy 



The differen t sources of error of our numerical method have been 
discussed in lValdettaro e t al. (20 071) and tested in a context sim- 
ilar to the present one in iLignieres et al.l (l2006h and lReese et al.l 
(l2006h . The numerical resolution has been chosen to ensure a 
sufficient accuracy for the computed frequencies. In the hori- 
zontal direction, the resolution is given by the truncation of the 
spherical harmonics expansion. The highest degree of the expan- 
sion is L = 2ne + \m\ and, for most of the calculations presented 
here, we used L - AQ + \m\, i.e. ng - 20 coupled spherical har- 
monics. In the pseudo-radial direction ^, the solution has been 
expanded over the set of Chebychev polynomials up to n^ - 96. 
Using higher resolutions (n^ - 96, L = 80 and n,. - 144, 
L - 80), we find that the relative agreement of the frequencies al- 
ways remains better than 5x10"**. It also does not affect the mode 
significantly as illustrated in Fig. |2] where the spectral expan- 
sion of the radial velocity component of the (f, m, n) = (1,0,14) 
mode at O = Q.1Q.K is displayed for the three different resolu- 
tions. Figure|2]also shows that a unique -or even a few- spherical 
harmonics would not properly describe such an eigenmode. 

2.4. Following modes with rotation 

We computed the frequencies of ^ = 1 to 3 modes in a nonrotat- 
ing polytrope. We recall that without rotation the system to solve 
becomes decoupled with respect to £, hence the modes are repre- 
sented with only one spherical harmonic. A reference frequency 
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Fig. 1. (top) Map in the meridional plane of the Brunt-Vaisala 
frequency No for the model with O = 0.70a:- No has been nor- 
malized by Q^ - sjGM/Rp where Rp is the polar radius. Dashed 
lines show the shape of the spheroidal grid we used, (bottom) 
Solid black lines show the profile of A^,, normalized by Q^ for 
the nonrotating star. Red lines show A^,, for the model Q = Q.1Q.K 
along the polar (solid line) and equatorial radius (dashes). A 
thick red tick on the x-axis indicates the polar radius Rp. 



set, w|,^', was computed from a 1-D poly tropic model with a ra- 
dial resolution n^ - 512. 

We then followed the variation in frequency of each mode of 
degree {q, azimuthal order mo, and radial order no by slowly in- 
creasing the rotation rate, step by step. The Arnoldi-Chebychev 
method requires an initial guess for the frequency, and returns 
the solutions that are the closest to this guess. The guess we pro- 
vide is extrapolated from the results at lower rotation rates: we 
compute from the three last computed points a quadratic extrap- 
olation at the desired rotation rate. For the first point (Q = 0), 
we use the frequency obtained in 1-D as a guess. For the sec- 
ond point, we extrapolate a guess with the asymptotic relation 
Wf„,,„„,„o ^ alfl,„^ + moQ/[A)(A) + 1)] (lLedouxlfT95l. 

Among the solutions found around the initial guess, we se- 
lect the correct one by following this strategy: 

1. For each calculated mode, we determine, from its spatial 
spectrum (like the ones shown Fig.|2]), the two dominant de- 
grees, £i and ti- 

2. We compare £i and £2 with the degree (q of the mode we are 
following. 
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Fig. 2. Spectrum C^ ^^^ ^ as a function of the degree ( of the ra- 
dial velocity of the mode (£,m,n) - (1,0, 14) at Q = OJD.K- 
^kmax^^ is the absolutc value of the highest Chebychev coeffi- 
cient in the decomposition on the spherical harmonics of degree 
( of the radial component of v. The spectrum is normalized to 
its maximum value. The different lines correspond to different 
spatial resolutions. 



3. We select the solutions such that {\ - {q; if none of the solu- 
tions verifies this criterion, we select the solutions such that 

ii - io- 

4. If more than one solution has been selected at this point, we 
consider the projection of the modes on the spherical har- 
monic Y"'° and compare it to the projection of the mode at a 
lower rotation rate. The solution that gives the highest corre- 
lation is finally selected. 

This method allows us to have a semi-automatic procedure 
that limits the need to manually tag the modes. Nevertheless, 
there are two limitations. The first one is inherent to the density 
of the g-mode spectr um. In deed, according to the asymptotic re- 
lation from .TassouH (Il980l) . on a given frequency interval, the 
number of modes of degree £ scales roughly as V^(^ + !)■ If our 
resolutions, both in latitude and radius, were infinite, it would 
almost be impossible to find the desired solution, beacause there 
would always be an infinite number of other solutions closer to 
the initial guess than the mode being searched. In practice, the 
spatial resolution is finite, and we noticed that, in most cases, the 
number of solutions in a small interval is limited enough to al- 
low us to find the desired solution. Another way of avoiding this 
difficulty is to add thermal dissipation that disperses the differ- 
ent solutions in the complex plane. This also proved successful 
in finding a specific solution. 

The second difficulty comes from the so-called avoided 
crossings. Two modes with the same m and the same parity can- 
not have the same frequency. This implies that the two curves 
associated to their evolution with Q cannot cross each other. 
Figure |3] illustrates this phenomenon with modes ^ = 1 and 
{ - 3: the frequencies get closer and closer, but since the curves 
cannot cross, the modes exchange their properties. During an 
avoided crossing, the two modes have the mixed properties of 
the two initial modes. With our mode-following method, when 
the coupling is strong and the avoided crossing takes long, the 
method can follow the wrong branch. For instance, in the case 
illustrated in Fig. |3] if the program follows the £ = I mode, it 
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Fig. 3. Evolution with the rotation rate of the frequency of the 
modes {I - \,m - \,n - \) and (/" = 3,m = l,n = 3), with a 
zoom on a region where an avoided crossing between these two 
modes occurs. At a given rotation rate, a diamond (respectively, a 
triangle) indicates the mode is dominated by the component ( - 
I it - 3). Crosses indicate that the modes are hardly discernible: 
they are both dominated by f = 1 and their structures are very 
similar 



continues sometimes on the f = 3 branch instead of jumping to 
the other branch. 



3. Perturbative coefficients 



and get 



^(,m,n + ^t,m,n^ + ' 



c?,t:y+o(^"') 



CL,„ + Ct,,„x + --- + d-^+0(x*-'). 



(15) 
(16) 



We note that Ci ,-, vanish for odd /. 

We compute y^^^ ^ and y^^ ^ on a grid of k points from Q = 
SO. to fc(5Q and use the Eqs. ( fTSl l and ( fTSI l to calculate the terms 
C^^jjj^ with the {k - l)-th-order interpolating polynomials. The 
determination of the coefficients Ci is then accurate to the 
(2k - l)-th-order in Q. In practice we use a typical step 5Q. - 



2 X 10"^ and A: = 4. In Eq. ([T4]i, we use 



:,(0) 



[wf_„,,„(Q 



0) + a)(-,„^„(il = 0)]/2, making it totally independent of the 1-D 
solutions. 

By explicitly expressing the dependence on m of the pertur- 
bative coefficients, Eq. (fTOb becomes 



- pert 



-,(0) 



miTl + m^TljD.'+Oid") (17) 

The form of the 1st or der comes from iLedouxl (Il95l|). the 2nd 
order from ISaiol (Il98lh . and the 3rd is derived from lSoufi et al.l 
(1199 8). We have verified that the derived coefficients fit these 
relations with a very good accuracy (see below) and list them in 
Tabled 

To know the coefficients for another normalization, for in- 
stance for d) - a>/Q.K, one can use the following development: 






(18) 



The approach used to determine the perturbative coeffi cients in The perfiirbed frequencies in this new normaUzation then ex- 
this paper is very close to the one of lReese et alj (l2006l) . press 



3.1. Determining perturbative coefficients 

In the perturbative approach, frequencies are developed as a 
function of the rotation rate, Q.. For instance, to the 3rd order, 
it reads 



a>' 



pert 



-m 



(0) 



„ = < + mC,„Q + {Sl^ + Aull + m'Sl)n'+ 



m{Tl^+m'Tl„)D.'+0{a''). (19) 



From our models we have computed A k 0.77164. 



'^On.n = ^S + Cl„,n^ + Cl„^„Cl^ + Cl^^„Cl' + 0{Cl\ (10) 



- pert 



where w~, ^ is the frequency for the non-rotating case, and C^ „ „ 
the perturbative coefficients. The bar denotes the normalization 
o) = ca/Q.^^ and Q - Q/il^. We normalize the frequencies by Q'^ 
since the polar radius is expected to be a slowly varying function 
of Q in real stars, as opposed to R^q. 

The coefficients C^^ can be numerically calculated from the 
complete computations since they are directly linked to the j-th 
derivative of the function <ijf,„,.„(Q) at Q = 0. However, to im- 
prove the accuracy, we use symmetry properties of the problem: 
changing Q in -Q., one easily shows that 



d ^(-lyci Vm. 



We define 



y£,m,n ~ 



ye,m,n 



2£1 



m > 0, 






m > 0, 



(11) 

(12) 
(13) 

(14) 



3.2. Coefficient accuracy and comparisons witli previous 
works 

The zeroth-order coefficients w^] were compared to the 1- 

D computations and we find agreement within 10"^. We 
also compared our results to previous frequency computa- 
tions of in a nonrotating polytrop ic model performed by 
IChristensen-Dalsgaard & Mullanl (IT994) with a totally diff'erent 
method. We renormalized their results for g modes (Table 4 of 
their paper) to their dynamical frequency Vg = 99.8557 yuHz (Eq. 
3.2 of their paper). The relative diff'erences with our results do 
not exceed 5 x 10"**. 

The choice for the step i5Q is important for the accuracy of 
the terms C^^^^^. Ideally, we should choose 6Q. as small as pos- 
sible, but when it is too small, the numerical noise, produced by 
the uncertainties on the computed (jL>e,m,n (Sect. 12.3b . drastically 
increases. We then chose the value of SO. to have the best trade- 
off. These uncertainties on Ci determinations were taken into 
account for the estimated accuracy of the coefficients C[^„, S'^^ 
and T'r . 

The Ist-order perturbative coefficients Cf,„ are expressed 
with integrals of the eigenmodes in the nonrotating model 
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Table 1. Perturbative coefficients (see development Eq. [TTl i for g modes with frequency u > 0.255, radial order n < 25, and { < 3 
in a poly tropic stellar model with an index fj. = 3. 



n 


<! 


Co. 




Si 




^ Cm 




{ = 1 












1 


1.5861677 


0.47187464 


0.0027 


-0.1215 


0.0742 


- 


2 


1.1338905 


0.46515116 


0.1946 


-0.0991 


0.1255 


- 


3 


0.8807569 


0.46565368 


0.3473 


-0.0779 


0.1769 


- 


4 


0.7195665 


0.46890412 


0.4801 


-0.0604 


0.2383 


- 


5 


0.6082150 


0.47259900 


0.6021 


-0.0458 


0.3131 


- 


6 


0.5267854 


0.47600569 


0.7176 


-0.0330 


0.4016 


- 


7 


0.4646791 


0.47895730 


0.8289 


-0.0215 


0.5037 


- 


8 


0.4157567 


0.48146486 


0.9374 


-0.0108 


0.6188 


- 


9 


0.3762235 


0.48358671 


1.0440 


-0.0007 


0.7469 


- 


10 


0.3436109 


0.48538634 


1.1490 


0.0089 


0.8877 


- 


11 


0.3162449 


0.48692017 


1.2530 


0.0182 


1.0411 


- 


12 


0.2929509 


0.48823507 


1.3562 


0.0272 


1.2069 


- 


13 


0.2728805 


0.48936906 


1.4586 


0.0360 


1.3852 


- 


14 


0.2554057 


0.49035278 


1.5606 


0.0445 


1.5758 


- 


1 


{ = 2 
2.2168837 


0.16413695 


-0.0603 


-0.1283 


0.1162 


0.0019 


2 


1.6817109 


0.13416727 


0.2040 


-0.1171 


0.1492 


-0.0186 


3 


1.3499152 


0.13379919 


0.3992 


-0.1267 


0.2263 


-0.0409 


4 


1.1271730 


0.13720207 


0.5686 


-0.1417 


0.3295 


-0.0666 


5 


0.9676634 


0.14084985 


0.7244 


-0.1587 


0.4555 


-0.0964 


6 


0.8478758 


0.14409412 


0.8717 


-0.1767 


0.6031 


-0.1307 


7 


0.7546269 


0.14685697 


1.0136 


-0.1952 


0.7716 


-0.1695 


8 


0.6799744 


0.14918725 


1.1517 


-0.2140 


0.9608 


-0.2128 


9 


0.6188542 


0.15115489 


1.2870 


-0.2330 


1.1704 


-0.2607 


10 


0.5678867 


0.15282450 


1.4201 


-0.2521 


1.4006 


-0.3132 


11 


0.5247312 


0.15424993 


1.5517 


-0.2713 


1.6500 


-0.3700 


12 


0.4877153 


0.15547464 


1.6819 


-0.2905 


1.9198 


-0.4315 


13 


0.4556123 


0.15653342 


1.8111 


-0.3098 


2.2095 


-0.4975 


14 


0.4275023 


0.15745412 


1.9395 


-0.3291 


2.5190 


-0.5679 


15 


0.4026818 


0.15825915 


2.0672 


-0.3484 


2.8482 


-0.6429 


16 


0.3806035 


0.15896664 


2.1943 


-0.3678 


3.1970 


-0.7223 


17 


0.3608352 


0.15959137 


2.3209 


-0.3872 


3.5655 


-0.8062 


18 


0.3430314 


0.16014547 


2.4470 


-0.4066 


3.9536 


-0.8945 


19 


0.3269120 


0.16063895 


2.5728 


-0.4259 


4.3613 


-0.9873 


20 


0.3122481 


0.16108017 


2.6983 


-0.4453 


4.7884 


-1.0846 


21 


0.2988504 


0.16147607 


2.8234 


-0.4648 


5.2351 


-1.1863 


22 


0.2865611 


0.16183254 


2.9484 


-0.4842 


5.7013 


-1.2924 


23 


0.2752477 


0.16215452 


3.0731 


-0.5036 


6.1869 


-1.4030 


24 


0.2647981 


0.16244623 


3.1976 


-0.5230 


6.6920 


-1.5180 


25 


0.2551166 


0.16271128 


3.3220 


-0.5424 


7.2164 


-1.6374 


1 


e = 3 

2.6013404 


0.06527813 


-0.1840 


-0.0702 


0.0898 


0.0117 


2 


2.0582624 


0.04834015 


0.0662 


-0.0551 


0.0478 


0.0014 


3 


1.6990205 


0.05125011 


0.2366 


-0.0576 


0.0593 


-0.0023 


4 


1.4466219 


0.05532661 


0.3818 


-0.0631 


0.0803 


-0.0054 


5 


1.2597371 


0.05898678 


0.5132 


-0.0695 


0.1071 


-0.0085 


6 


1.1157943 


0.06207639 


0.6359 


-0.0764 


0.1389 


-0.0120 


7 


1.0015072 


0.06465872 


0.7527 


-0.0836 


0.1753 


-0.0158 


8 


0.9085566 


0.06682318 


0.8653 


-0.0909 


0.2162 


-0.0200 


9 


0.8314693 


0.06864897 


0.9747 


-0.0983 


0.2616 


-0.0246 


10 


0.7664974 


0.07020016 


1.0817 


-0.1058 


0.3113 


-0.0296 


11 


0.7109883 


0.07152736 


1.1869 


-0.1133 


0.3652 


-0.0350 


12 


0.6630115 


0.07267052 


1.2906 


-0.1208 


0.4235 


-0.0409 


13 


0.6211284 


0.07366127 


1.3931 


-0.1284 


0.4859 


-0.0471 


14 


0.5842454 


0.07452488 


1.4946 


-0.1360 


0.5526 


-0.0537 


15 


0.5515158 


0.07528169 


1.5953 


-0.1437 


0.6234 


-0.0608 


16 


0.5222741 


0.07594820 


1.6953 


-0.1513 


0.6985 


-0.0683 


17 


0.4959899 


0.07653788 


1.7947 


-0.1590 


0.7776 


-0.0761 


18 


0.4722351 


0.07706183 


1.8936 


-0.1666 


0.8609 


-0.0844 


19 


0.4506608 


0.07752926 


1.9920 


-0.1743 


0.9483 


-0.0931 


20 


0.4309792 


0.07794783 


2.0901 


-0.1820 


1.0399 


-0.1022 


21 


0.4129512 


0.07832398 


2.1878 


-0.1897 


1.1356 


-0.1117 


22 


0.3963764 


0.07866313 


2.2853 


-0.1974 


1.2353 


-0.1216 


23 


0.3810855 


0.07896988 


2.3825 


-0.2051 


1.3392 


-0.1320 


24 


0.3669346 


0.07924815 


2.4794 


-0.2128 


1.4472 


-0.1427 


25 


0.3538006 


0.07950128 


2.5762 


-0.2205 


1.5592 


-0.1538 
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(lLedouxl|l95lh . We then computed these terms with our 1-D 
eigensolutions and compared them to C^(„,„/m. The results are 
consistent within 10"'^. 

An explicit computation of 2nd- and 3rd-order coefficients 
requires calculating the 1st- and 2nd-order corrections of the 
eigenfunctions, which is not so straightforward. It is the rea- 
son we performed a direct numerical computation of these co- 
efficients. The numerical errors we estimated for S', and Tl 

are generally around 10"^ and always less than lO"'*. We then 
checked the co nsistency of o ur computations with the 2nd-order 
calculations of I Said ( Il981h for g modes with n - 1 to 3. In 
this work, all frequencies were normalized by the dynamical 
frequency Q.j/ of the nonrotating polytrope. By noticing that 

Q/Q^" = Q H- AQ.^ + 0{Cl^) with A * 0.18391, and using the 
relation iT% . we were able to compare these results with ours. 
We get a good qualitative agreement with absolute differences 
better than 10"^, which is reasonable relative to the lower accu- 
racy of Saio's computations. It gives an interesting consistency 
check for our calculations. Overall, the perturbative coefficients 
listed in Table [Tlhave been determined with high accuracy. 



4. Domains of validity of perturbative approaches 

From the previously computed coefficients, we calculated mode 
frequencies with the 1st to 3rd-order perturbative approxima- 
tions for rotation rates ranging from Q = to 0.7Q/(- and com- 
pared them to complete computations. Figure |4] illustrates such 
a comparison by showing the evolution of the frequencies of 
the seven m components of an /' = 3 mode together with their 
2nd-order perturbative approximation. We clearly observe that 
the agreement between both approaches at low rotation progres- 
sively disappears as the rotation increases. 
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Fig. 4. Evolution with the rotation rate of the frequencies of the 
components of the ({ - 3, n - 16) multiplet obtained with a 
complete computation (solid line) and with the 2nd-order per- 
turbative approximation (dashed lines). 

To define the domains of validity of perturbative approaches, 
we fix the maximal departure 6a> allowed between the perturbed 
frequencies wj?^^^'^^ and the "exact" ones Wf_„,„. For each mode 
and each approximation order, we define the domain of validity 
[0,Q,,], such that |w^'"„(n) - Wf,,„,„(Q)| < (5a> VQ < Q,. The 
precision of the observed frequencies 6v can be related to the 



6a> — 2n 6v 



GM' 



(20) 



From ffiis expression, we see that, for a fixed precision 6v, 
the normalized error 5cj depends on the dynamical frequency 
Vg = {GMIR^)'- l{2n) of the star considered. We thus display 
the domains of validity of the perturbative approximations for 
two types of stars with different dynamical frequencies, a typi- 
cal y Dor star, and a typical B star The y Dor star is such that 
M - 1.55Mo, R - 1.6Rq, i.e. Vg - 61 iiHz, while the B star has 
M = 4-Mq, R = IRq, i.e. v^ = 1 1 jjHz. For the frequency preci- 
sion, we chose 6v -Q. lyuHz, which corresponds to the resolution 
of an oscillation spectrum after a hundred days. This is the typi- 
cal accuracy for a CoRoT long run. Accordingly, the normalized 
error Sd) is equal to 1.6 x 10"^ for the y Dor star and 9.3 x 10"^ 
for the B star. It must be noted that the domains of validity will 
not be affected by the numerical errors of wj?"^ ^ (cf. Sect. 13.21 ) 
and a)(j„^„ (cf. Sect. 12.31 ) as the values of doj remain at least an 
order of magnitude higher 

We have determined the domains of validity of 1st-, 2nd-, 
and 3rd-order methods for low-degree { < 3 modes. Specifically, 
we considered { = 1 modes with n = 1 to 14, and { = 2 and 3 
low-order (n = 1 to 5), and high-order (n = 16 to 20) modes. 
The domains of validity are shown in Fig. |5]for both types of 
stars. Overall, the domains of validity extend to higher rotation 
rates for B stars than for y Dor stars. This is simply due to the 
increase in the normalized tolerance 6a). Besides, we observe 
distinct behaviors in the high- and low-frequency ranges. 

In the high-frequency range, 2nd-order perturbative meth- 
ods give satisfactory results up to ~100 kms"' for y Dor stars 
and up to ~150 km s ' for B stars. The 3rd-order terms improve 
the results and increase the domains of validity by a few tens of 
km s"'. These results are to be contrasted with those found for p 
modes where the domains of validity are restricted to lower ro- 
tation rates. F or 6 Scuti stars, wh ich have similar stellar param- 
eters to y Dor. lReese et al.l (12006*) find ~50-70 km s ' as a limit 
for perturbative methods. In addition, the 3rd-order terms do not 
improve the perturbative approximation in this case, as p modes 
are weakly sensitive to the Coriolis force. The rather good per- 
formance of perturbative methods at describing high-frequency 
g modes indicates in particular that the 2nd-order term gives a 
reasonable description of the centrifugal distortion. This might 
be surprising considering the significant distortion of the stellar 
surface (Req/Rp - 1.08 at Q - OA^k)- Actually, the energy 
of g modes is concentrated in the inner part of the star where 
the deviations from sphericity remain small (as shown in Fig.[T} 
top). As a result, g modes "detect" a much weaker distortion 
that is then amenable to a perturbative description. A particu- 
lar feature that induces a strong deviation from the perturbative 
method concerns mixed pressure-gravity modes that arise as a 
consequence of the centrifugal modification of the stellar struc- 
ture. For example, we found that, above a certain rotation rate, 
the ^ = 3,n = 1 mode becomes a mixed mode with a p-mode 
character in the outer low-latitude region associated with a drop 
in the Brunt-Vaisala frequency No (see Fig. [TJ. 

The domains of validity of perturbative methods are strongly 
reduced in the low-frequency range. For yDor stars, 2nd-order 
perturbative methods are only valid below ~50 kms"'. Indeed, 
a striking feature of Fig. |5]is that perturbative methods fail to 
recover the correct frequencies in the inertial regime w < 2Q 
(delimited by a magenta curve). In particular, we observe that, 
although increasing the tolerance Sai between the left (y Dor) 
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Fig. 5. Evolution of the frequencies of /" - 1,2,3 modes (top to bottom). Frequencies are computed in the corotating frame. 
Perturbative approximations have been tested for a typical -y Dor (left panels) and for a B star (right panels). Green/red/blue parts 
of curves indicate that lst/2nd/3rd order is sufficient to reproduce complete calculations within an error 6v - Q.l //Hz. Error bars 
on the righthand side of each panel show 6v and lOx 6v. Magenta lines indicate a> = 2Q. For each plot, the bottom x-axis and left 
y-axis show dimensionless units, whereas the top x-axis and right y-axis show physical units. 



and right (B star) panels subtantially extends the domains of va- in the mode cavity that are not taken into account by the per- 
lidity in the a> > 2Q regime, very little improvement is observed turbative method. Indeed, we observed that modes in the iner- 
in the o) < 2Q regime. tial regime do not explore the polar region and that the angular 

size of this forbidden region increases with 2Q/w. This is illus- 
In the following, we argue that the failure of the perturbative tj-^ted in Fig. i for a particular mode. Such a drastic change in 
method in the subinertial regime [0, 2Q] is related to changes 
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Fig. 6. (Left) Meridional distribution of kinetic energy jPi,v^ of 
the g mode {( - 3,m - -\,n - 16) in a nonrotating star To 
enhance the contrast, it is scaled by r^. (Right) The same for 
Q - 0.70a:- Magenta lines indicate the critical surface F = 0. 

the shape of the resonant cavity has a direct impact on the asso- 
ciated mode frequency. As perturbative methods totally ignore 
this effect, they cannot provide an accurate approximation of the 
frequencies in this regime. 

This interpretation is supported by t he analytical expres- 
sion of the forbidden region determined bv lDintrans & RieutordI 
( 12000 ) for gravito-inertial modes. Indeed for frequencies co < 
2£l, modes are mixed gravity-inertial modes, since the Coriolis 
force becomes a restoring force. In the context of their spherical 
model, and within the anelastic approximation and the Cowling 
approximation, they have shown that gravito-inertial waves with 
a frequency at only propagate in the region where 



rV[A?2 + (2Q)2. 



w"] - (2Q.Nozf > 0. 



(21) 



This implies that, when w < 20, a critical latitude 0^ - 
arcsin[(y/(20)] appears above which waves cannot propagate. 
Even though this expression does not strictly apply to our non- 
spherical geometry, we have overplotted the critical surfaces 
r = with the energy distributions of our eigenmodes (see Fig. |6] 
for an illustration). For the nonrotating case, there is only a small 
circle close to the center, corresponding to the classical turning 
point 0) = No- For the mode with w < 20, the polar forbidden 
region delineated by F = agrees pretty well with the energy 
distribution of our complete computation. 

5. Conclusion 

In the present work, we have computed accurate frequencies 
for g modes in poly tropic models of uniformly spinning stars. 
We started from high- and low-frequency, low-degree (^ < 3) g 
modes of a nonrotating star and followed them up to fi = Q.1Q.K- 
This allowed us to provide a table of numerically-computed per- 
turbative coefficients up to the 3rd order for a polytropic stellar 
structure (with index ^ - 3). This table can serve as a reference 
for testing the implementation of perturbative methods. We then 
determined the domains of validity of perturbative approxima- 
tions. For the high-frequency (low-order) modes, 2nd-order per- 
turbative methods correctly describe modes up to ~ 100 kms ' 
for y Dor stars and up to ~ 150 kms ' for B stars. The do- 
mains of validity can be extended by a few tens of km s"' with 
3rd-order terms. However, the domains of validity shrink at low 
frequency. In particular, perturbative methods fail in the inertial 
domain o) < 2Q because of a modification in the shape of the 
resonant cavity. 

In a next step, we plan to compare our complete com- 
putations with the so-called traditional approximation, which 
is also extensively used to determine g-mode frequencies (e.g 



iBerthomieu et al.lll978HLee & Saidll997h . We will also analyze 
how rotation affects the regularities of the spectrum - such as 
the period spacing - and compare it to the predictions of the per- 
turbative and traditional methods. In the present study, we have 
focused on low-degree modes, but a more complete exploration 
clearly needs to be performed. I n particular, we might look for 
the singular modes predicted bv iDintrans & RieutordI (12000*). It 
requires to take care of dissipative processes, which play an im- 
portant role in this case. 
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